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q^' Abstract 

rn | The corrections for two pion correlations due to electromagnetic final state inter- 

actions at high secondary multiplicities are investigated. It is shown that these result 
in a noticeable deviation from the standard Gamov factor. This conclusion changes 
drastically in a model of the pion system with expansion. The critical parameter which 
determines the size of these effects is found to be the ratio of the relative velocity of 
detected pions to the velocity of the pair center-of-mass (in the fireball rest frame). In 
particular, when this parameter is much less than unity the pion pair escapes the initial 
high density region promptly and the distortion of the mutual Coulomb potential is 
weak. 
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The fundamental observable for intensity interferometry in hadron physics is the relative 
momentum spectrum between two identical particles. For two like-charged pions, the final- 
state Coulomb interaction modifications to this spectrum result in a correction which has 
typically been considered to be tractable and relatively accurate. This expectation is based 
on the significantly different length scales between the strong (oc l/m n ) and the Coulomb 
(oc l/m n a) interactions in the problem ^ (here m n is a pion mass and a is a fine structure 
constant). In this case the the correction may be treated with the Schrodinger equation, 
resulting in the well-known Gamov factor G(pi,p2). The nominal quantity expressing the 
correlation function in terms of experimental distributions [^J 
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n(n \ < n > d 3 Pl d 3 p 2 m 
C(Pl ' P2)= <n(n-l)>rfVrfV' (1) 

d 3 pi d 3 p 2 

where < n > is the particle multiplicity, d 3 a/d 3 p and d 6 a / d 3 p\d 3 p2 are the single-pion and 
two-pion cross sections, is then given (due to the factorization of the corresponding matrix 
element) in terms of the product of a Gamow factor multiplying the model correlations j3J 

C(Pi, P2) = C(pi, p 2 )C model (p 1 , p 2 ). (2) 

The standard derivation of the Gamov factor [Q, obtained from the solution ip(r) of the 
nonrelativistic Schrodinger equation with the Coulomb potential, leads to 



g{q) =1 m r= ( 3) 



where 



V = Vo= -g-, (4) 

where Q is the relative momentum between the particles. It should be mentioned that, 
in terms of the pion momentum k in the pair center-of-mass, the relative pion momentum 
Q = 2k coincides with invariant relative momentum Qi nv = [{k\ + /E2) 2 — 4m 2 ] 1 / 2 where 
k\ and k2 are pion four-momenta in an arbitrary frame. In order to better illustrate the 
more complicated cases discussed below, it is worthwhile to reconsider the calculation of this 
phenomenon as a quantum mechanical tunneling process (see for instance 0). 

In the quasi-classical approximation the probability for a pion of mass starting from 
the point r2 to reach the distance r\ under the barrier is related to the quantity 

1 /" r 2 

r)pen{ri,r 2 ) = - / q(r)dr, (5) 



7T Jri 

with 



q(r) = [2m 7V (v(r)-E° km )} 1/ \ (6) 

being the so-called "subbarrier quasi-momentum" . The equation V(r 2 ) = E® in determines 
the classical turning point. 
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For the pure Coulomb barrier Vcw = ct/2r the calculation of Eq. (jS) gives 



Vpenin, r 2 ) = - [I(r 2 ) - J(ri)] , (7) 

7T 



with 



/(r) = rq(r) — rjo arcsin 



■y 2 kin 



V Co ui{r) 



where the factor of 1/2 in the potential accounts for the usage of the physical pion mass, 
not the reduced one, in Eq. @. For penetration to r\ = 0, corresponding to production of 
the two pions at the same space-time point, one obtains 

V P en(ri = 0, r 2 ) = = 77o, (9) 

identical with (|). 

Turning to the high multiplicity case, the relation between the two-particle electromag- 
netic potential and the local charge density is given by 

V 2 0(r) = -4vre(n (+) - n ( ' y ) , (10) 

where e is the elementary charge and and where the density of charged mesons n^' is then 
related to that of neutral mesons via a Boltzmann factor: 

n (±) = n (°) exp , (11) 

Here is the density of 7r°-mesons at the freeze-out temperature T, which coincides with 
the equilibrium density of charged pions in the absence of the Coulomb interaction (we 
consider symmetrical nuclear matter). When e<p <C T the Eq.(|TT|) can be rewritten as 

n& = n<°) (l^) , (12) 

(this requires that the pions are not closer than ~ 10~ 2 fm to one another at T « 200 MeV) 
so that 

Attp 2 

V 2 0(r) = _(2n<°>Mr) , (13) 

where e 2 = a. 

We can write down immediately the well-known solution of the Eq.flTB]) as the screened 
Coulomb potential 

7r± ( r ) = ±e (14) 

r 

where 



a ■ \ — , (15) 
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with potential energy U nn = a exp(— r / R scr ) / r for the like-sign pions. 

To evaluate the correction to the Gamow factor, we will use this screened Coulomb 
potential to re-calculate the penetration parameter given by Eq.([5]). This in turn requires 
an estimate of the pion density and freeze-out temperature. We will calculate this for the 
extreme case, i.e., just after freeze-out, when pions occupy the volume |7[] (we take the same 
temperature Tf = 190 MeV) 

V f = irRpTf sinh , (16) 

where for Pb-Pb collisions Rf = Rpt> ~ 7.1 fm, Tf = 10 fm and Ay « 6. For multiplicities 
we assume |§ : LHC: N n = 8000, SPS: N w = 800, which gives for screening radii LHC: R scr = 
7.9 fm, SPS: R scr = 25 fm . The latter looks rather controversial at first sight (though 
for a pure Coulomb interaction R scr = oo). Indeed, using Eq. ([UJ) for the unit interval of 
rapidity (one might consider that more natural, in a sense) we obtain much smaller value of 
corresponding screening radius. None the less we cite this quantity here as an upper bound 
in existing experimental conditions. 

The results of these calculations together with the standard Gamov factor are plotted 
in Fig. 1. Also shown there is a curve obtained by direct substitution of the experimental 
values obtained by NA44 |§ into Eq. (|1§), with T f = 187 MeV, r f = R L « 6.0 fm, 
Rt — 6 fm, Ay = 3 and dN/dy = 40. These parameters produce an even smaller screening 
radius (R scr — 19.3 fm) than our nominal SPS value (and hence a larger correction to the 
standard Gamov factor), which follows from the smaller values of R scr ~ 4^jTf/n n at the 
NA44 experimental conditions. Based on these consideration alone, it would appear that a 
substantial correction to the Gamow factor is already required by the existing experimental 
data. This is particularly true for kaon interferometry since the pion medium screens the 
K-K Coulomb final state interaction as well, and the length scale for the K-K Gamov factor 
is larger than that for pions by the ratio of their masses. 

However, it is important to note that correction factors presented in Fig. 1 are an upper 
bound, in that they do not incorporate the expansion of the pion source after freeze-out. We 
next turn to a more realistic calculation explicitly incorporating expansion. It is intuitively 
clear that the correction factor for the fast pairs will approach the standard Gamov factor 
and for the very slow ones the estimate obtained should be valid. 

We parametrize the spherical expansion of the pion source 

«(fl)=»/§F. ( 17 ) 

in terms of the parameters rif, the freeze-out pion density; the corresponding radius Rf 
and the freeze-out temperature Tf. In this model the spatial volume of the expanding pion 
system in the solid angle Q increases as Q ■ R 2 , where R is the distance from the center of 
the fireball. Then, the corresponding potential is the solution of the Maxwell equation 

(V 2 - df) 0(r) = ^n^(r) , (18) 

where we put rS^ w n n /3 as before. Now introducing the distance r between two detected 
particles we have 
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R Vrrn. ' t 



Vrel ■ t 



(19) 



where v cm is the velocity of the two-particle center-of-mass in the fireball rest frame and 
v re i is the relative velocity of the particles (v re i = Q/m, t = is fixed on the freeze-out 
hyper-surface). Then 



R 



r . 



Vrel 



so that Eq.(|T%D may be rewritten as 



(V 2 - v 2 re tf) 0(r) 



c\Q) 



with c 2 {Q) defined as 



c\Q) 



87re 2 RfU^ 



rel 



T vl 



(20) 



(21) 
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Fixing the screening radius at the freeze-out temperature and density as 

R f = 



3T ( 



/ 



we have 



and Eq.(IT^) can be rewritten as 



c(Q) 



d 



57r an f 



Rf Vrel 
Rscr Vcm 



dr 

This equation has the following solution 



r 2 (l-v 2 rel )^ + 2r^-c\Q) 



dr 



. 



eRo 



a-l 



(23) 



(24) 



(25) 



(26) 



where Rq is a scale parameter to fix the dimensions of the potential 0. The exponent a is 
a solution of the quadratic equation resulting from the substitution of the ansatz (f2^) into 
Eq . (|2l~D and takes the form (assuming proper asymptotic behaviour of the potential </>) 



1 

a = - 
2 



l + vf 



rel 



1-v 2 



+ 



rel 



l+V, 



rel 



Vrel, 



4c 2 (Q) 
1-v? 



rel 



(27) 



For small pion relative velocity v re i (corresponding to Q ~ — 30 MeV) the above equation 
reduces to the simple expression 



l + 4c 2 (Q) . 



(2* 
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Both are shown for several sets of parameters in Fig. 2. In the interval of interest "a" , and 
hence deviations from a pure Coulomb field, increases with increasing relative velocity. This 
intriguing behavior is directly related to the "Hubble-like" expansion implied by Eq. (|20|). 
and becomes quite understandable if we remember that the pion density decreases (hence 
R scr — > oo) with R increasing. The expansion thus results in modifications to the Coulomb 
potential that are of power-law, not exponential form, in contrast to that static result given 
by Eq. ([Hp. We can treat the corrected potential obtained in Eq.(PB[) in terms of an effective 
charge distribution 



'eff 



Rq 

r 



(29) 



which we are going to average over using the following procedure. If we confine out attention 
to v re i <C 1, Eq.(|21|) can be rewritten in a Poisson-like form (V 2 — K 2 )<f)(r) = 0, with 
k = c(Q)/r, which is equivalent to a r-dependent screening radius, i.e., 



c(Q) 



(30) 



Since, as shown in Fig. 2, the deviation of the potential ( [26D from the pure Coulomb form 
in the region of small relative pion momenta Q < 30 MeV is small, we first ignore the r 
dependence of k to obtain the solution for the electromagnetic tttc potential in the form of 
Eq. (|iHD, then substitute the r-dependent n into Eq. (fTJ) to find 



IK 



ae 



-c(Q) 



(31) 



Instead, the numerator on r.h.s. 
squared (we are now considering 
b, hence the extra factor of charge 



Thus, there is no an exponential dependence of U nn on r. 
of Eq.(|3lD represents the averaged charge distribution ( [29^ 
the potential energy U nn rather than the electric potential 
leading the a). 

It is clear from Eq. fl2"4|) that when the ratio v re i/v cm <C 1 (Rf and R* scr are of the 
same order for high multiplicities) the renormalized constant aexp[— c(Q)] is close to the 
bare value of a. Moreover, the same qualitative result comes from the r-behaviour of the 
screening radius ( |3"0"D when it approaches the asymptotic value R scr = oo (Coulomb law) with 
increasing r. The quantity c(Q) increases with relative pion momentum leading to larger 
deviations from the Coulomb potential and agreement with the features of the potential (p6|) 
(see discussion after Eq. (|28D). 

Based on these considerations, solving the Schrodinger equation for the screened potential 
gives results similar to Eqs. (131) and (f§), but with the renormalized a, so that 



am 



-exp[—c(Q)] = i] ■ exp 



Q 



and taking into account the Eq. (|15D we have 



Vrel R 



f 



Vcm Rscr 



(32) 



Rf_ 

11.SCT 



1 2a 



R f T f 



(33) 
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This result shows explicitly that the large modifications to the Gamow function from to 
screening are weakened in the expansion model due to the small value of v re i/v cm . Specific 
examples of this are shown in Figs. 3 and 4. 

Since for much of the parameter space c(Q) <C 1, it is useful to further approximate 
exp[— c(Q)]. Again referring to the NA44 data ||, and considering the direction transverse 
to the collision axis where v cm ~ p?/ % for the set of parameters 



p T > 150 MeV/c, N n < 200, R f > 4 fm, T f > 180 MeV, Q < 40 MeV/c, (34) 



one has exp(— c(Q)) ~ 1 — c(Q), and the screening of the Coulomb interaction reduces to 
shift rj 



7] ~ r] — Ar] 



with 



Ar] = a — — ■ * 2a — — 
Pt V T f R f 



(35) 



(36) 



being independent of relative pion momentum Q. (Clearly, for high multiplicities and for 
small pt this approximation will be violated.) 
For the case of cylindrical geometry we have 



and for the ratio 



Rf_ 
p/ 

-Ttscr 



rrRpTf sinh \Ay 



(37) 



a 



(38) 



and finally 



rj = r/o ■ exp 



2QrriT 



aN v 



PtTU^ \ 3TjTf sinh \Ay 
Again, for the set of parameters where one might expect this geometry to apply, 



(39) 



p T > 150 MeV/c, N v < 2000, 77 > 10 fm, Ay > 3, T f > 180 MeV, Q < 40 MeV/c . 

(40) 

we obtain the shift At] as 



TTIt 

Arj = 2a— 



aN v 



Pt \ 3TfTj sinh ^Ay 
which is also independent of relative pion momentum. 



(41) 
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Fig. 3 shows the role of screening for the correction factor G scr behaviour evaluated at 
several sets of parameters together with the standard Gamov factor G . Their ratio G /G scr 
plotted in Figs. 4 and 5 show that correlators measured in the transverse direction develop 
a well-pronounced dependence on the transverse momentum of the pair in the region of 
Pt < 150 MeV/c JO]] where the deviation from standard Gamov factor behaviour increases. 

Our considerations show that for future LHC and RHIC experiments the screening radius 
of Coulomb interaction at the freeze-out density and temperature could be comparable with 
the source size and therefore the factorization of Eq. (0) |4| is no longer valid. However, 
this conclusion changes drastically with the inclusion (switching on) of expansion for the 
pion system. We would like to emphasize that this main result could be model independent. 
The detailed evaluation of the Gamow correction in pion interferometry analysis at very 
high multiplicities of secondary particles in the picture of an expanding fireball reveals an 
important regulating parameter what is the ratio of relative velocity of the detected pions 
and their center-of-mass velocity in the rest frame of a fireball v re i/v cm . 
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Figure 1. Extreme correction factors for screened Coulomb potential versus relative 
pion momentum. Bottom curve is the standard Gamov factor. 
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Figure 2a. Exponent a in the deformed Coulomb law </>(r) = ei?Q _1 /r a , 
which was calculated with taking into account the dependence of <f> on time 
(see Eqs.(28) and (30)) for the following set of parameters: 
1 - R f = 6 fm, N v = 1200, p T = 450 MeV/c, T f = 190 MeV; 
■2 (i ////. .V^ 120, p T = 150 MeV/c } T } = 187 \l< \ : 

3 - R f = 4 fm, N v = 120, p T = 450 MeV/c, Tj = 190 MeV. 




Figure 2b. Exponent a in the deformed Coulomb law </>(r) = ei?Q _1 /r a 
which was calculated without taking into account the dependence of <f> on 
time (see Eq.(31)) for the same as in Fig. 2a set of parameters. 
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Figure 3. Correction factor for screened Coulomb potential versus relative 
pion momentum (expansion scenario). G scr was evaluated for the set of 
parameters: 

1 - R f = 6 fm, N v = 1200, p T = 450 MeV/c, T f = 190 MeV; 

2 - R } = 4 fm, N v = 120, p T = 450 MeV/c } T } = 190 MeV; 

3 - standard Gamov factor. 




Figure 4. Ratio of the standard Gamov factor and correction factor for 
screened Coulomb potential versus relative pion momentum (expansion sce- 
nario). G scr was evaluated for the set of parameters: 

1 - R f = 4 fm, N v = 120, p T = 450 MeV/c, T f = 190 MeV; 

2 - LHC conditions r s = 10 fm, N v = 8000, p T = 150 MeV/c } Ay = 6, 
T f = 190 MeV; 

3 - R f = 6 fm, N v = 1200, p T = 450 MeV/c, T } = 190 MeV. 




Figure 5. Ratio of the standard Gamov factor and correction factor for 
screened Coulomb potential versus relative pion momentum (expansion sce- 
nario). G scr was evaluated for the set of parameters: 

1 - R f = 6 fm, N v = 120, p T = 150 MeV/c, T f = 187 MeV; 

2 R } 1 ////. AT„ 120, p T = 50 MeV/c, T f = 187 \l< \ : 
:i R } :i ////. A.. 30, ,>, = 20 MeV/c, T } = 187 MeV. 



